Strong Field Ionization Rate for Arbitrary Laser Frequencies 
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A simple, analytical, nonrelativistic ionization rate formula for atoms and positive ions in intense 
ultraviolet and x-ray electromagnetic fields is derived. The rate is valid at arbitrary values of 
the Keldysh parameter and confirmed by results from ab initio numerical solutions of the single 
active electron, time-dependent Schrodinger equation. The proposed rate is particularly relevant for 
experiments employing the new free electron laser (FEL) sources. 
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In the year 2000, when a free electron laser (FEL) came 
into operation at DESY [ij, a new era of strong field 
laser physics has begun: the interaction of intense, ultra- 
violet and x-ray laser pulses with matter. Presently, in- 
tensities ~ lO^^W/cm^ are available at photon energies 
Huj ~ 100 eV 0. One of the most fundamental quan- 
tum processes in intense laser-matter interaction is the 
nonlinear photoeffect where several photons are absorbed 
by the emitted electron. Ionization of atomic targets by 
intense pulses from FELs has been studied in several ex- 
periments (3, S, II i, H, the most recent of which [2] 
reports ionic charge states up to Xe^^+ at 13.3 nm and a 
laser intensity of about 10^^ W/cm^. 

The photoeffect becomes highly nonlinear and requires 
a nonperturbative description if Kq = I /huj 3> 1, with / 
the ionization potential. For intense, low-frcqucncy fields 
where ionization proceeds via tunneling through the bar- 
rier formed by the binding potential plus the quasistatic 
laser field such nonperturbative descriptions have been 
developed decades ago. According to the standard termi- 
nology, the tunneling regime is determined by the small- 
ness of the Keldysh parameter 01, 7 = V^w/^o ^ 1, 
where is the electric field amplitude of the laser wave 
(we use atomic units unless noted otherwise). In this 
limit an analytical expression for the total ionization 
rate of atoms and positive ions was found for the first 
time in Ref. Q. Later this result was extensively gen- 
eralized, including relativistic ionization and the ioniza- 
tion of molecules or other spatially extended systems 
(see the review |9|, also Refs. [l^ [llj, and references 
therein). Compared to tunneling out of a short-range 
potential, the long-range Coulomb interaction between 
the outgoing electron and the atomic core of charge 
Z > \ strongly increases the tunneling probability be- 
cause the Coulomb attraction suppresses the potential 
barrier through which the photoelectron tunnels. Ex- 
perimentally, the Coulomb-induced enhancement of the 
tunneling ionization rate is well established since the 80s 
[13] ■ Nowadays the Coulomb-corrected tunneling rates 
(widely known as "ADK rates" [H, [li]) are commonly 



used for the calibration of laser pulse intensities. Much 
less attention has been paid to the theory of the opposite 
multiphoton limit 7 3> 1. For strong fields this regime 
can only be reached using high frequencies, which is the 
reason why it became experimentally accessible only with 
the invention of the high-power FELs. In the experi- 
ments [1, H, 0, @ typical values of the Keldysh parameter 
were 7 ~ 30 100. Tunneling rates are completely ir- 
relevant in this domain: at an intensity ~ lO^^W/cm^ 
(as used in Ref. for instance, sequential tunneling 
ionization would predict ions up to Xe^^-Xe^"*" only in- 
stead of Xe^^+, and the probability to ionize Xe^°+ via 
tunneling is ~ 10^^^ for a lOfs laser pulse. 

In this Letter an analytical expression for the total 
(energy and angle-integrated) ionization rate of atoms 
and positively charged ions in the field of an intense, 
high-frequency laser pulse, i.e., for 7 ^ 1, is derived, 
including the Coulomb-correction, which gives an up to 
nine orders of magnitude enhancement. Our rate yields 
the correct tunneling limit so that it is actually valid at 
arbitrary frequencies and arbitrary values of the Keldysh 
parameter. The only limitations are that (i) our rate 
is restricted to the nonrelativistic regime, and (ii) that 
only one electron is considered active. We prove the high 
accuracy of our result by comparing the predicted rates 
with ah initio numerical results. 

The ionization rate w of an atom (ion) with the elec- 
tron in a state of ionization potential / can be represented 
as the product of the rate wsr for a system bound by 
a short-range potential with the same ionization poten- 
tial and the Coulomb correction Q, which accounts for 
the asymptotic Coulomb interaction Uc ~ —Z/r, r 3> 
I/V2I between the outgoing electron and the core [l5| . 
The theoretical approach for the derivation of Q based 
on the method of classical, complex trajectories (imag- 
inary time method [lB|) was introduced in Q for the 
tunneling regime 7 <C 1. Readers interested in the de- 
tails of the complex trajectory method are referred to 
Refs. 0, [13, 0, m [iSl where also several applications 
are presented. In the semiclassical limit Kq 3> 1 the ion- 
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ization probability is determined by the imaginary part of 
the reduced classical action evaluated along a trajectory 
in complex time: 



exp{-2ImVK}, = / {C - I) dt-wv 
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Here, C is the Lagrangian, r(t), v(t) are the electron tra- 
jectory and velocity, respectively, and is the complex 
start time determined by the condition v^(ts) = —21. A 
trajectory r(i) satisfies Newton's equation of motion. If 
only the laser force is accounted for in this equation [we 
denote a corresponding Coulomb- free trajectory as ro(0]i 
Eq. ([1]) gives the probability w^k for the ionization from 
a short-range potential. This is the well-known output of 
Keldysh theory 0, @] or the strong field approximation 
2lj . A perturbative account of the Coulomb force yields 



the Coulomb correction 



Q = cxp{-2ImWc}, 



(2) 



where F = Eq/{2I)'^/'^ is the reduced electric field and 
n* = Z/V2T is the effective principle quantum number 
of the bound state. The ionization rate of an atom (ion) 
then is ^ 



W = Q ■ WSR, 
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Here, r{t) ~ ro{t)+ri{t) is the Coulomb-corrected trajec- 
tory, S£ = VqVi +vJ/2 — f (t) -ri is the respective correc- 
tion to the Lagrangian, 5(v-r) = Vo-ri-|-Vi-ro-l-Vi-ri, and 
£{t) is the electric field of the laser pulse. To eliminate 
the divergences in ([3]) a matching with the atomic wave 
function is required. This matching technique and other 
computational details are described in Refs. 17, 1^. 

In order to find the correction to the total ionization 
rate it is sufficient to consider the most probable electron 
trajectory. In a linearly polarized field £{t) = £0 cos{ujt) 
the photoelectron momentum vanishes at infinity (i.e., 
at the detector) for this most probable trajectory. In 
the limit 7 3> 1 the respective Coulomb-free trajectory 
xo(t), the correction xi{t) (both one-dimensional, along 
the polarization direction) and the complex start time ts 
have the form [17] Xo{t) — 6(1 ~ u), u — e'^^'^", < r = 
—iujt < tq, 

xi{T) = b{i^/2iLT + ^[T^/2-\n{l-u) + L2(u)]}. (4) 



Here b = v2//a;, tq = —iujts ~ ln(27), and L2{u) is the 
Euler dilogarithm [2^. The parameter fi ~ Zlo j {2Tf'l'^ 
determines the relative contribution of the Coulomb field 
to the trajectory (at 7 ^ 1) so that a perturbative ac- 
count formally requires, in particular, <C 1 (see the 
discussion of the applicabily conditions in [l^l)- Note 
that because of the term ~ ^JJi in (|4]) we keep the seem- 
ingly higher-order term vi • ri in (S(v - 0, appearing in 
([3|) (see the corresponding remark in |20[). After some 
cumbersome but straightforward algebra one obtains for 
the Coulomb-correction ^ 
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e = 2.718. 
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ci = arcsinh 7 — 7/^/1 4- 7^, /? = 2'y/yjl + 7^, T{x) is 
the Dawson integral |23j[, and nth — Kq[1 + 1/(27^)] is 
the ionization threshold (in units of the photon energy). 
For the asymptotic coefficient C of the bound state wave 
function we use Hartree's approximation [S^l, = 



-)2n*-2 



/{n*\f (for s-states). 



Expression ^ is our main result. Although derived 
for 7 ^ 1, the Coulomb correction ([5]) obeys the right 
tunneling limit. Thus for 7 ^ 1 Eq. ([6]) gives the well- 
known expression for the tunneling ionization rate of an 
atomic s-state [1, . In the intermediate domain 7 ~ 1 
our result can be considered as an interpolation. In the 
limit 7^1 the correction is intensity-independent and 
numerically large, Q ~ (2eiiro)^" ^ 1- For noninteger 
nth the Dawson integral is of the order of unity, the rate 
([7]) simplifies for 7^1, and one obtains for ([6]) 
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where iVm = [Kq] + 1 is the minimum number of photons 
required for ionization. The fact that the rate ^ is pro- 
portional to the intensity to the power of is, of course, 
not surprising in the multiphoton regime where an atom 
absorbs the minimum possible number of light quanta 
0, @| . The new achievement is the derivation of the co- 
efficient A(KQ,n*) with, as we show below, high quanti- 
tative accuracy. Note that the coefficient A{Ko,n*) is a 
product of big numbers and thus very sensitive to small 
variations. The Coulomb correction (O affects the value 
of A significantly: for, e.g., the parameters Kq ~ 5 10, 
n* ~ 1-^3, typical for strong-field ionization of atoms 
and positive ions, Q varies between 10^ and 10^*^. The 
second factor in Eq.®, which distinguishes our Coulomb 
correction from the tunneling limit, alters the ionization 
rate also significantly for 7 ^ 1. Note that in Ref. 15 [ 
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FIG. 1: Ionization rate (in atomic units) vs laser intensity 
according ((6| (solid), ([7|) (dashed), and from the respective 
tunneling formula (dotted line). Rates obtained numerically 
using the Floquet method and from the solution of the TDSE 
are shown by dots and triangles, respectively. Upper panel: 
H(ls) in a field corresponding to the second harmonic of a 
Ti:Sa leiser (A = 400 nm). Insert: same but for A — 422 nm. 
Lower panel: Xe^'^^ (4po-electron with / — 434 eV, n* = 3.19) 
in a A = 13.3-nm x-ray field. 



the Coulomb-free probability wsn was calculated for ar- 
bitrary 7 but the Coulomb correction was taken in the 
tunneling limit Q = {2/F)'^"\ 

As examples, the rates (jH]) and ([7]) are plotted as a 
function of the laser intensity in Fig. [T] for (i) the ion- 
ization of atomic hydrogen by the second harmonic of 
a Ti:Sa laser and (ii) for the ionization of Xe^^+ by a 
13.3-nm x-ray field. The Coulomb-parameter fi is 0.11 
for H and 0.34 for Xe so that the theory is applicable 
in both cases. The correction ^ is about 6 x 10^ and 
9 X 10^, respectively. The corresponding tunneling rate 
is included in both panels, showing that it is completely 
irrelevant for 7 3> 1 but merges with ^ when approach- 
ing the tunneling limit 7 <§C 1. In order to check the 
analytical expressions we compare their predictions with 
exact numerical results obtained by two different meth- 
ods: the Floquet method yields the exact ionization rate 
for an infinite pulse as twice the (absolute value of the) 
imaginary part of the Floquet-eigenenergy. The solution 
of the time-dependent Schrodinger equation (TDSE) in 
real time in general depends on the pulse shape but with 



increasing pulse duration also yields converged results 
for the rate. We used the publicly available codes STR- 
FLO and Qprop In the real-time TDSE calcula- 
tions the norm Af{t) inside a sufficiently large test-sphere 
around the ion is calculated. The ionization probability 
1 — A/'(oo) divided by the pulse duration r is then taken 
as the ionization rate. Alternatively, one may determine 
the rate from the slope of A/'(t). The 4po state of Xe^''"'" 
is constructed as an eigenstate of the effective potential 
U{r) = -Z/r - (54 - Z)e-'"7r, where Z = 18 and the 
parameter k = 7.93 is adjusted to reproduce the ioniza- 
tion potential / = 434 eV. 

The comparisons of our analytical rate with the nu- 
merical results show, in general, a good, quantitative 
agreement even for the case of xenon where the pa- 
rameter fi = 0.34 is not so small and the bound state 
is a p-state while the analytical result (l5])-(l7]) is de- 
rived for s-states. Deviations of our analytical rate from 
the numerical results can be attributed to resonances. 
This is shown in the upper panel of Fig. [T] where a 4- 
photon resonance with a Stark-shifted Rydberg level is 
met. This resonant enhancement of the rate at intensi- 
ties ~ 2 3 X 10^^ W/cm^ is more pronounced in the 
Floquet results because this method assumes an ideally 
monochromatic pulse. In the real-time TDSE solution a 
10-cycle flat-top pulse with a two cycles up- and down- 
ramp is used, which has a finite spectral width and thus 
smears-out resonances. As it is seen from the insert the 
analytical result agrees very well with both sets of nu- 
merical data when the resonance is avoided by a small 
variation of the wavelength. 

Besides our results in this regime we are only aware of 
the ones obtained by Santra and Greene for the ioniza- 
tion of Xe at Huj — 12.7 eV and an intensity 10^"^ W/cm^ 

f, which were the parameters used in the experiment 
where charge states up to Xe^+ were observed. In 
Ref. 27[ the multiphoton cross-sections — wn/]^ 
were found numerically, with TV having the same meaning 
as A^in in Eq. ([5]) and j — cSq/Sttlo the photon flux. The 
same cross-sections can be easily calculated from Eq. ([5]) . 
A comparison between our results and the cross-sections 
of Ref. 27[ show, in general, reasonable agreement. 
For example, we find from Eq. ([5]) 0-4 (Xe'^"'") ~ 2.8 x 
lO-^^cm^s^ and CT5(Xe'*+) ~ 5.3 x 10-i5°cmi0s4, while 
Santra and Greene give tT4(Xe^"'") ~ 1.8 x 10~^^^cm^s^ 



and CT5(Xe 



1.1 x 10" 



respectively. Our 



analytical cross-sections are systematically lower, a fact 
we attribute to the method used for the extraction of the 
rate from the numerical data. More precisely, a complex 
absorbing potential for distances r > Rq — A a.u. was 
used in 27| so that population of Rydberg states may 
have contributed to ionization. 

A simple analytical expression for the multiphoton ion- 
ization rate for Rydberg states of atoms and ions was 
derived long ago by Berson [2^ under the condition 
n* > while our case is the opposite, N^ji ~ Kq > n*. 
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As a consequence, the result of Berson yields a different 
expression for the coefficient A{Ko,n*) in Eq. ([8]) and 
does not match the tunneling limit. 

In conclusion, we presented an analytical formula for 
the ionization rate of atoms and ions which is valid for ar- 
bitrary values of the Keldysh parameter. The high accu- 
racy of our rate was demonstrated by comparisons with 
exact numerical calculations. In the present form, our 
theory discards resonances and describes ionization from 
s-states in a linearly polarized field. Generalization to 
the case of arbitrary polarization and arbitrary angular 
momentum of the initial state, and the inclusion of rela- 
tivistic corrections are straightforward. 

Our theory is based on the single active electron ap- 
proximation. Thus in complex atoms it should de- 
scribe the sequential ionization channel when electrons 
are emitted one after another. It is by now well es- 
tablished, however, that the electron-electron interac- 
tion may change the ionization rate dramatically via 
correlated, nonsequential mechanisms. In low-frequency 
fields this nonsequential ionization is essentially under- 
stood within the recollision or 'atomic antenna' picture 
(2^ . [sol . 31 1 . In the high-frequency domain nonsequen- 



tial mechanisms have been studied both in experiments 
and theoretically for cases where onl y a few laser pho- 



tons are involved (see, e.g. |33j, [SJ] and references 

therein). In the multiquantum domain we consider here 
multielectron effects remain to be explored. The theo- 
retical work by Santra and Greene, e.g., indicates that 
the cross-sections change up to one order of magnitude 
if many-electron effects are taken into account 27|. The 
experiment reported in Ref. ^2] indicates that for highly 
charged xenon ions the double-logarithmic slope of the 
intensity-dependent yield is not equal to but rather 
remains constant with increasing charge state. These 
facts are hints that multielectron mechanisms may also 
play an important role in multiquantum ionization at 
high frequencies. However, even if this is the case our 
rate is useful since any measured deviation from the pre- 
dictions of Eqs. ([5])~([7]) which do not sensitively depend 
on the wavelength (and therefore must be nonresonant in 
nature) can then be attributed to electron-electron cor- 
relation. 
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